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Abstract. We report on the nonlinear turbulent processes as- 
sociated with electromagnetic waves in plasmas. We focus 
on low-frequency (in comparison with the electron gyrofre- 
quency) nonlinearly interacting electron whistlers and non- 
linearly interacting Hall-magnetohydrodynamic (H-MHD) 
fluctuations in a magnetized plasma. Nonlinear whistler 
mode turbulence study in a magnetized plasma involves in- 
compressible electrons and immobile ions. Two-dimensional 
turbulent interactions and subsequent energy cascades are 
critically influenced by the electron whisters that behave 
distinctly for scales smaller and larger than the electron 
skin depth. It is found that in whistler mode turbulence 
there results a dual cascade primarily due to the forward 
spectral migration of energy that coexists with a backward 
spectral transfer of mean squared magnetic potential. Fi- 
nally, inclusion of the ion dynamics, resulting from a two 
fluid description of the H-MHD plasma, leads to several 
interesting results that are typically observed in the solar 
wind plasma. Particularly in the solar wind, the high-time- 
resolution databases identify a spectral break at the end of 
the MHD inertial range spectrum that corresponds to a high- 
frequency regime. In the latter, turbulent cascades cannot be 
explained by the usual MHD model and a finite frequency 
effect (in comparison with the ion gyrofrequency) arising 
from the ion inertia is essentially included to discern the dy- 
namics of the smaller length scales (in comparison with the 
ion skin depth). This leads to a nonlinear H-MHD model, 
which is presented in this paper. With the help of our 3D 
H-MHD code, we find that the characteristic turbulent in- 
teractions in the high-frequency regime evolve typically on 
kinetic-Alfven time-scales. The turbulent fluctuation asso- 
ciated with kinetic-Alfven interactions are compressive and 
anisotropic and possess equipartition of the kinetic and mag- 



netic energies. 
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1 Introduction 

Many laboratory and space plasmas contain multi-scale elec- 
tromagnetic fluctuations. The latter include the electron 
whistlers and H-MHD fluctuations in a uniform magneto- 
plasma. While the wavelengths of the electron whistlers 
could be comparable with the electron skin depth, the H- 
MHD fluctuations could have differential wavelengths when 
compared with the ion skin depth and the ion-sound gyrora- 
dius. In whistler mode turbulence, the electrons are incom- 
pressible and the ions form the neutralizing background in 
the plasma. On the other hand, in the H-MHD turbulence the 
compressibility of the electrons and ions cannot be ignored. 

In uniform magnetized plasmas, one can have different 
types of electromagnetic waves. The latter include the 
low-frequency (in comparison with the electron gyrofre- 
quency) electron whistlers as well as linearly coupled fast 
and slow magnetosonic and kinetic-Alfven waves. The elec- 
tron whistlers are described b y the electron magn etohydro- 
dynamics (EMHD) equations Kingsep et al ( 1990), and the 
rol e of whistlers in the EMHD turbulence is an unre s olved is- 
sue IShaikh et al I ri2000albh: IShaikh & Zank I d2003l): IShaikh 
"IT 



d2004l) : IShaikh & Zank I d2005h : IShaikh I (120081 12009b . ~No 
tably, E MHD in two dimensio ns exhibits dual cascade phe- 
nomena Biska mp et al. ( 1996|) simil a r to th ose of 2D Navier- 
Stokes dKraichnanl 1 19651; iFrischl 1 1995b and MHD dy- 
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namic s dFvfe & Mon tgomery T~ 1977 : iBiskamp & Bremer 
19941) . However, because of the complexity of nonlinear 



interactions in EMHD turbulence, mode coupling interac- 
tions are relatively more complex to comprehend. Like 2D 
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Navier-Stokes dKraichnanl 1 19651: iFrischL 1 19951) an d MHD 



(iFyfe & Montgomery I 



1977 



Biskamp & Bremer L 11994 



systems, 2D EM HD also exhibits dual cascade processes 



(Bis kamp et al. 



1996). The dual cascade phenomena, i.e. 
both the forward and inverse cascades, is however diffi- 
cult to resolve numerically in a single energy spectrum be- 
cause of a lack of spectral resolution in the inertial ranges. 
Therefore there exists considerably less numerical work on 
this issue as compared to various two-dimensional (decay- 
ing) turbulence systems. Numerical simulations of inverse 
cascade phenomenon, within the context of magnetohy- 
drodynamics (MHD) turbulence, have previously been per- 
formed at a modest resolu tion of up to 1024 2 Fourier modes 
dBiskamp & Bremer L[l994h . 

The dynamics of linearly coupled magnetic field-aligned 
Alfven wave and obliquely propagating fast and slow 
magnetosonic waves and kinetic Alfven waves are governed 
by the H-MHD equations. The solar wind plasma wave 
spectra corresponding to whistlers or high-frequency (in 
comparison with the ion gyrofrequency) kinetic-Alfven 
wave possess a spec tral break, the origin of which is not 
yet fu lly understood Denska tet al dl983h: Goldstein et al 



also re ported near the Earth's magnetopause IStenberg et al 
(2007) and the Cluster spacecraft encountered them during 
the process of magnetic reconn ection in the Earth's mag- 
netotail region IWei et al I ( 12007b . Upstream of collisionless 
shock, whistler wave s are found to play a c rucial role in heat- 
ing the plasma ions IScholer & Burgess I d2007l) . Their ex- 
citation and propagation are not only limited to the Earth's 
nearby ionosphere, but they are also found to be excited 
near the ionosphere of other planets such as in the radia- 
tion belts of Jupiter and Saturn iBespalov I (|2006). Whistlers 
are believed to be a promising candidate in transporting 
fields and cu r rents i n plasma opening switch (POS) devices 
iMason et al I d 19931) . which operate on fast electron time 
scales. Similarly, electron whistlers have been known to 
drive the phenomen on of magnetic field line r econnection 



Bula nov et al I d 19921) in astrophysical plasmas Zho u et al 
dl996l) . Whistlers have also been investigated in several labo 



ratory e xperime nts Stenzel ( 1997 l|l975l) : Istenzel & Urrutia 



>ry e 

(1990); iStenzel et al I (11993.199401) . where they have been 



jl994: iLeamonetall dl998l): 



Goldstein et al 



dl996l) : 



Bale et al. ( 2005 ): Servidio et al. I d2008h . There are, by rec ent three-dimensional simulations [Eliasson & Shukla 



found to exhibit a variety of interesting features, such as 
an anisotropic propagation of the phase front, strong dis- 
persion characteristics, interaction with plasma particles, 
etc. A few experimental features have also been confirmed 



however, some pr e liminary res ults available in the literature 
Shaikh & Shuklal d2008ll2009l) . 

In this paper, we present a survey of the spectral properties 
of the low-frequency multi-scale electromagnetic turbulence 
in plasmas. In section 2, we begin by describing the whistler 
wave model based on two dimensional electron magnetohy- 
drodynamic equations. Linear properties of whistler waves 
are described in this section. Section 3 contains nonlin- 
ear simulation results describing the wave spectra in small 
and large scale (in comparison with the electron skin depth) 
regimes. Section 4 describes the Hall- MHD counter part. 
The summary and conclusions are contained in section 5. 



2 The Whistler Mode Turbulence 

In an incompressible electron-MHD (E-MHD) plasma, col- 
lective oscillation of electrons exhibit electron whistlers. 
The electron Whistlers are widely observed in many space 
and laboratory plasmas. For instance, they are believed to 
be generated in the Earth's ionospheric region by lightning 
discharges and proceed in the d irection of Earth's dipole 
magnetic field (Hel liwell . 119651) . They have also been re- 
cently detected in the Earth's radia t ion b elt by the STEREO 
S/WAVES instrument ICattell et all d2008l) . Moreover, there 
are observations of Venus' ionosphere that reveal strong, 
circularly polarized, electromagnetic waves with frequen- 
cies near 100 Hz. The waves appear as bursts of radia- 
tion lasting 0.25 to 0.5s, and have the expected properties 
of whistler-mode si gnals genera t ed by lightning discharges 
in Venus' clouds iRussell et al I ((2007). These waves are 



(2007), where it has been reported that the polarity and 
the amplitude of the toroidal magnetic field, in agreement 
with the laboratory experiments, determine the propaga- 
tion direction and speed of the whistlers. These are only 
a few examples amongst a large body of work devoted 
to the study of the whistlers. Despite the large amount 
of effort gone into understanding the existence and the 
propagation of whistlers, their linear dynamics is still de- 
bated, specially in the context of complex nonlinear pro- 
cesses. For example, the role of whistlers in the high- fre- 
quency turbulence and anisotrop i c spectr al c ascades has been 
debate d recently [sh aikh et al (2000a b); Shai kh &Zank ' 



d2003l) : IShaikhl d2004l) : IShaikh & Zank I d2005l) . In the fol- 
lowing, we present a summary of the spectra of the whistler 
turbulence. 

The electron whistler mode dynamics is essentially gov- 
erned by the EMHD equations. The E MHD phenomena typ i- 
cally occur on the electron time scale IKingsep et all dl990)l. 



while t he ions do not participate in the whistler mode dy- 
namics IShukla I dl978l) l. Thus, the basic frequency regimes 
involved are uj C i <C u> -C 0J C e (where w C i,uj ce are, respec- 
tively, the ion and electron gyrofrequencies, and ui is the 
characteristic whistler frequency), and the length scales are 
c/uj p i < t < c/uj pe , where u> p i,ujp e are the ion and elec- 
tron plasma frequencies, respectively. In whistlers, electrons 
carry currents, since immobile ions merely provide a neu- 
tralizing background to a quasi-neutral EMHD plasma. The 
electron momentum equation is 



m e n—V e — — enE V e x B VP - /im e nV e , (1) 

dt c 
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Fig. 1. (left) Driven 2D whistler wave turbulence with the forcing band identified by the kf modes. The spectral resolution is 5120 2 . 
Notice that no spectral break exists across the forcing modes, which agrees with the Kolmogorov theory. Across the kf modes, the forward 
cascade of energy (k > kf) and the inverse cascade of mean magnetic potential (k < kf) yield identical spectral indices and are close 
to Ek ~ k~ 5 ^ 3 . We plot the compensated energy spectrum k 5 ^ 3 Et ~ e 2 ^ 3 ~ Const to show inertial range modes along a straight 
line that is indicative of constant energy dissipation rates in steady state driven dissipative turbulence, (right) The inverse and forward 
cascade regimes are identified in 2D driven whistler wave turbulence. The compensated energy spectrum in the inverse cascade regime 
k 7/3 E inv ~ e 2/3 fe 7/3 fc" 1 ~ e 2/3 k 4/3 , whereas the forward cascade is given by k 7 ^ 3 Ef w d ~ Const. The approximate error in the spectral 
slope is about ±0.072. The spectral break, separating a forward and an inverse cascade regimes, corresponds to a hump above the forcing 
modes kf that results from a direct excitation of the Fourier modes in driven whistler wave turbulence. 
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c~dt' 



„ „ 4vr 18E 
V xB = — J + - — , 

c c at 



On 

a* 



V • (nV e ) = 0, 



(2) 
(3) 
(4) 



where 



±_d_ 

dt~ dt + e 



The remaining equations are B = V x A, J = — enV e , V • 
B = 0. Here m e ,n, V e are the electron mass, density and 
fluid velocity respectively. E, B respectively represent elec- 
tric and magnetic fields and cj>, A are electrostatic and elec- 
tromagnetic potentials. The remaining variables and con- 
stants are, the pressure P, the collisional dissipation /i, the 
current due to electrons flow J, and the velocity of light c. 
The displacement current in Ampere's law Eq. © is ig- 
nored, and the density is considered as constant throughout 
the analysis. The electron continuity equation can, therefore, 
be represented by a divergence-less electron fluid velocity 
V • V e = 0. The electron fluid velocity can then be associ- 
ated with the rotational magnetic field through 

-V x B. 



47me 



By taking the curl of Eq. (Q]) and, after slight rearrange- 
ment of the terms, we obtain a generalized electron momen- 
tum equation in the following form. 



dP 

--V e x(VxP) 



where 



m P V» 



eA 

c 



and £ = -m e V e ■ V e 



(5) 



P 



The curl of Eq. © leads to a three-dimensional equation 
of EMHD describing the evolution of the whistler wave mag- 
netic field, 



d_ 

dt 



ft. 



V e • Vfl B - fl B ■ VV e = /<V 2 B. 



(6) 



where = B — dgV 2 B,<i e = c/oj pe , the electron skin 
depth, is an intrinsic length-scale in the EMHD plasma. The 
three-dimensional EMHD equations can be transformed into 
two dimensions by regarding variation in the z-direction as 
ignorable i.e. d/dz — 0, and separating the total magnetic 
field B into two scalar variables, such that B = z x + bz. 
Here ip and b respectively present perpendicular and parallel 
components of the wave magnetic field. The corresponding 
equations of these components can be written in a normalized 
form as follows, 



d_ 

dt 



z x V6 • Vfi 



B 4 b 

ay 



0. 



(7) 
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—to B -d 2 J x V6-VV 2 6+z x Vi/'-VVV+Sot^V 2 ^ = 0. 
at ay 



where toy = ip - d 2 V 2 ^,ri B = b - d 2 e V 2 b The length 
and time scales are normalized respectively by d e and u> ce , 
whereas magnetic field is normalized by a typical mean B . 
The linearization of Eqs. (|7]l & ® about a constant magnetic 
field Bq yields the dispersion relation for the whistlers, the 
normal mode of oscillation in the EMHD frequency regime, 
and is given by 



ky k 

i + dik 2 



where u Co = eB /mc and k 2 = k 2 + k 2 . From the set of the 
EMHD Eqs. (0 & (0, there exists an intrinsic length scale 
d e , the electron inertial skin depth, which divides the entire 
spectrum into two regions; namely short scale (kd e > 1) 
and long scale (kd e < 1) regimes. In the regime kd e < 1, 
the linear frequency of whistlers is ujf. ~ k y k and the waves 
are dispersive. Conversely, dispersion is weak in the other 
regime kd e > 1 since Uk ~ k y /k and hence the whistler 
wave packets interact more like the eddies of hydrodynami- 
cal fluids. 



3 Energy spectra in whistler turbulence 

By virtue of d e , there exists two inertial ranges that corre- 
spond to smaller kd e > 1 and larger kd e < 1 length-scales 
in the EMHD turbulence. Correspondingly, forward cascade 
turbulent spectra in these regimes exhibit fc~ 5 / 3 and fc~ 7 / 3 
respectively as shown in Fig. (Q~|i. The whistler spectrum 
£,-7/3 j s p ro( j ucec i essentially by fluctuations in electron fluid 
while ions are at rest. The observed spectra in our simula- 
tions depicted as in Fig. (fl]i can be explained on the basis of 
Kolmogorov phenomenology as follows. 

The regime kd e > 1 in EMHD turbulence corresponds 
essentially to a hydrodynamic regime because the energy 
spectrum is dominated by the shorter length-scale turbulent 
eddies that give rise to a characteristic spectrum of an in- 
compressible hydrodynamic fluid. The group velocity of 
whistlers, in this regime, is small and hence it is expected 
that the whistler effect cannot be present. For kd e > 1, the 
first term in the magnetic potential can be dropped and thus 

A ~ (V 2 ^) 2 . 

Similarly the dominating terms in the energy are 

E ~ (V6) 2 + (W 2 tp) 2 . 
This implies 

E ~ (V6) 2 ~ (V 2 i/0 2 . 

Thus for kd e > 1, E ~ A. In this case both magnetic poten- 
tial and the energy cascades turn out to be identical leading to 



(8) the same spectral index, i.e. A; -5 / 3 , in both direct and inverse 
cascade regimes. This is shown in Fig. (Q} (left panel). 

In EMHD the eddy velocity (v) in the x — y plane is char- 
acterized by z x V6. Thus the typical velocity with a scale 
size t can be represented by v ~ be,/ 1. The eddy scrambling 
time is then given by 



v bt 

In the limit of kd e < 1, the mean square magnetic potential 
is, 

A~i, 2 ~b 2 l 2 . 

The second similarity follows from the assumption of an 
equipartition of energy in the axial and poloidal components 
of magnetic field. This assumption has been invoked in some 
earlier works too, and would be exact if the EMHD tur- 
bulence comprised of randomly interacting whistler waves 
only. To determine the scaling behaviour of the spectrum in 
the magnetic potential cascade regime (i.e. for k < kf, kf 
representing the forcing scales), the rate of A transfer is given 
by the relation, 

F - — - h 3 
T 

and the mean square magnetic potential per unit wave num- 
ber is 



A,, 



2 y.3 



bit 



The locality of the spectral cascade in the wavenumber space 
yields Ak = e a k^. On equating the powers of bi and £, 
we obtain a = 2/3 and (3 = —3. In kd e < 1 regime, the 
expression for the energy can be approximated as E ~ b 2 ~ 
(V<0) 2 . This further implies that E = A/£ 2 and 



E, 



A similar analysis in the other regime of A: > kf where en- 
ergy cascade is local in wave number space leads to an energy 
spectrum with the index of —7/3 such that 



E k 



k~ 7 '\ 



The compensated spectra of Fig. (Q~|i (right panel) are consis- 
tent with above description. 

In the following section, we shall include the ion dynamics 
that correspond to turbulence scales in which kdi > 1, where 
di is ion skin depth. This regime of the magnetofluid dynam- 
ics is referred to as the Hall MHD, and it stems from the 
combined motion of the electron and ion fluids amidst pre- 
dominant density fluctuations that evolves on Kinetic-Alfven 
wave (KAW) time (and length) scales. 
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4 The High-Frequency Kinetic-Alfven Inertial Regime Hall-MHD equations, given by 



The high-frequency kinetic-Alfven regime of the plasma re- 
sults in the solar wind turbulence when plasma wave fluctu- 
ations comprise a multitude of time and space scales that are 
shorter than the ion gyroperiod (27r/f2 c i) and comparable (or 
longer) to (than) the ion thermal gyroradius/ion-sound gyro- 
radius or shorter than the ion skin depth, respectively. These 
electromagnetic fluctuations cannot be d escribed by the usual 



magnetohydrodynamic (MHD) model (Denskat et al , 1983 



Goldstein etail ll994t ILeamon et all 1 19981: iGoldstein et al 



19961: Shaikh & Shukla U2008[l2009h which involve the char 



acteristic time scale much longer than the ion gyroperiod. 
Typical ly, reconnection of ma gnetic field lines near the solar 
corona (Bhattacharjee , 2004), Earth's magnetosphere and in 
many laboratory plasma devices (IBhattacharjee et al I 2001) 
are other avenues where collisionless or kinetic processes 
call for inclusion of the finite Larmor radius and finite fre- 
quency effects in the standard MHD model. In particular, in 
the context of the solar wind plasma, higher time resolution 
databases identify a spectral break near the end of the MHD 
inertial range spectrum that corresponds to a high-frequency 
(> f2 c ,) regime where turbulent cascades are not explain- 
able by Alfvenic cascades. This refers to a secondary iner- 
tial range where turbulent cascades follow a fc~ 7 / 3 (where 
k is a typical wavenumber) spectrum in which the charac- 
teristic fluctuations evolve typically on kineti c -Alfve n time 
scales dHowes et al. L 120081: iMatfhaeus et al. L E008). The 
onset of the second or the kinetic-Alfven inertial range 
still eludes our und e rstanding of solar wind fluctu ations 
( ILeamon et all Il998t Ishaikh & Shukla L l2008i l2009h . The 
mechanism leading to the spectral break has been thought 
to be either mediated by the kinetic- Alfven waves (KAWs), 
or damping of ion-cyclotron waves, or dispersive processes, 
or other implicit nonlinear interactions associated with ki- 
netic effects in solar wind plasmas jPenskat et al L 1 1 983l 
Goldstein et al L 1 1994 ILeamon et al L I 19981: IGoldstein et al L 
19961: IShaikh & Shukla L 120081 l2Q09h . Motivated by these 
issues, we have developed three dimensional, time depen- 
dent, compressible, non-adiabatic, driven and fully paral- 
lelized Hall-magnetohydrodyn amic (H-MHD) simulations 
(Shaikh & Shukla , 2008, 20095 to investigate turbulent spec- 
tral cascades in collisionless space plasmas. 

In the high-frequency regime, u> > C j, the inertialess 
electrons contribute to the electric field which is dominated 
essentially by the Hall term corresponding to J x B force. 
The latter, upon substituting in the ion momentum equation, 
modifies the ion momentum, the magnetic field and total en- 
ergy in a manner to introduce a high-frequency (uj > il C i) 
and small scale (k±pL > 1, where is ion thermal gyrora- 
dius) plasma motions. The characteristic length scales (fc _1 ) 
associated with the plasma motions are smaller than the ion 
gyroradius (pl)- The quasi-neutral solar wind plasma den- 
sity (p), velocity (U), magnetic field (B) and total pressure 
{P = P e + Pi) fluctuations can then be cast into a set of the 



| + V.(pU) = 0, 



d 



di 



dB 

~dt ~ 

de 

di + ' 

where 



e = 



1 



p — + U • V U = -VP + -J x B 



V x U x B 



J x B 



r ? V 2 B, 



(9) 
(10) 

(11) 



-pU 2 U 



p 

(7^1) 



— - — —p\J + — E x B ) =0(12) 

7 — lp 4-7T 



8tt 



is the total energy of the plasma that contains both the 
electron and ion motions. All the dynamical variables are 
functions of three space and a time, i.e. (x,y,z,t), co- 
ordinates. Equations Q to (fl"2l are normalized by typ- 
ical length Iq and time to = £o/vo scales in our sim- 
ulations, where vo — Bq/IAttpq) 1 / 2 is Alfven velocity 
such that V = l V,d/dt = t d/dt,tj = U/w ,B = 
B/wo(47rpo) 1 ^ 2 , P — P/po v q,P = p/po- The parameters 
p and r\ correspond respectively to ion-electron viscous drag 
term and magnetic field diffusivity. While the viscous drag 
modifies the dissipation in plasma momentum in a nonlin- 
ear manner, the magnetic diffusion damps the small scale 
magnetic field fluctuations linearly. The magnetic field is 
measured in the unit of Alfven velocity. The dimension- 
less parameter in magnetic field Eq. (fTTl i i.e. ion skin depth 
di = di/£o,di = C/uj p i is associated with the Hall term. 
This means the ion inertial scale length (di) is a natural or an 
intrinsic length scale present in the Hall MHD model which 
accounts for finite Larmour radius effects corresponding to 
high frequency oscillations in kdi > 1 regime. Clearly, 
the Hall force dominates the magnetoplasma dynamics when 
l/p(J x B) > U x B term in Eq. (fTTb which in turn 
introduces time scales corresponding to the high frequency 
plasma fluctuations in kdi > 1 regime. Furthermore, our 
model includes a full energy equation [Eq. (fT2l) l unlike an 
adiabatic relation between the pressure and density. The use 
of energy equation enables us to study a self-consistent evo- 
lution of turbulent heating resulting from nonlinear energy 
cascades in the solar wind plasma. 

In the linear regime, the Hall-MHD equations (without the 
energy equation) admit a dispersion relation, 

{u 2 - k 2 z V 2 )D m (io,k) = u 2 {u 2 - k 2 V 2 )k 2 k 2 V*/Lo 2 cl , 

which exhibits coupling among the fast and slow magne- 
tosonic waves given by the solutions of 



-D m (w,k) 



u 2 k 2 {vi + v;) + Kk 2 viv 2 = o, 
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where co is the frequency, k(= + k z z is the wave 
vector, the subscripts _L and z represent the components 
across and along the external magnetic fields Bqz, and Va 
and V s are the Alfven and effective sound speeds, respec- 
tively. The warm HMHD plasma thus supports a great va- 
riety of waves (e.g. magnetic field-aligned non-dispersive 
Alfven waves, obliquely propagating fast and slow magne- 
tohydrodynamic waves, obliquely propagating dispersive ki- 
netic Alfven waves and whistlers) having different wave- 
lengths (comparable to the ion skin depth associated with the 
Hall drift, the ion sound gyroradius associated with the elec- 
tron pressure and the perpendicular ion inertia/ion polariza- 
tion drift). Interestingly, in the w < Vl C i regime , which is pre 



domin antly known to be an Alfvenic regime, iHowes et al. 
(2008) noted the possibility of occurance of highly oliquely 
propagating KAWs (with uo <C f2 c i) thereby questioning the 
role of damping of ion cyclotron waves in the onset of spec- 
tral breakpoint. While our simulations described in Fig. © 
may contain the KAWs corresponding to u> <C f2 c i> their spe- 
cific contribution in the spectral cascades depicted in Fig. © 
and Fig. (0 is not identifiable. We therefore leave this study 
to a future investigation. 

The nonlinear spectral cascades in the KAW regime lead to 
a secondary inertial range in the vicinity of modes kdi ~ 1 
in that the solar wind magnetic, velocity and field fluctua- 
tions follow turbulent spectra close to fc~ 7 / 3 . This is shown 
in Fig. (0. The characteristic turbulent spectra in the KAW 
regime is steeper than that of MHD inertial range which fol- 
lows a k~ a where a is 5/3 or 3/2. The onset of the sec- 
ondary inertial range is constantly debated because of the 
presence of multiple processes in the KAW regime that in- 
clude, for instance, the dispersion, damping of ion cyclotron 
waves, turbulent dissipation or due to fast or slow magne- 
tosonic perturbations and etc. In the context of our simula- 
tions, the observed fc -7 / 3 in the solar wind turbulent plasma 
can be understood from the energy cascades affected by the 
Hall forces. The latter are one of the potential candidates that 
may be responsible for a fc~ 7 / 3 spectrum in the KAW inertial 
range regime. 

To understand the observed k~ 7 ^ 3 spectrum in the KAW 
inertial range regime in solar wind plasma and in our sim- 
ulations, we in voke Kolmog o rov and Kraichnan like p he- 
nomenologies (iKolmogorovL 1 1 95 lb iKraichnan I 1 1 9651) in 
Hall MHD plasma. The time scale associated with Hall 
MHD is th smaller than that of MHD tmhd- The non- 
linear cascades in the MHD turbulence are governed typi- 
cally by tmhd ~ (fcufc) -1 , where is the velocity field 
in the /c-space. By contrast, the spectral transfer of the 
turbulent energy in the solar wind (Hall MHD) plasma is 
th ~ (fc 2 £?fc) -1 . The energy transfer rates in the KAW 
regime are therefore e ~ u\/th. On substituting the turbu- 
lent equipartition relation between the the velocity and mag- 
netic fields B\ ~ uf,, the energy transfer rates turn out to be 
e ~ k 2 B\. The use of the turbulent equipartition is justi- 
fied from the observations and it is also seen in our 3D sim- 




Fig. 2. Inertial range turbulent spectra of magnetic and velocity field 
fluctuations are shown. These fluctuations closely follow a k~ 7 ^ 3 
spectrum in kinetic Alfven regime where kdi > 1. The simulation 
parameters are: spectral resolution is 128 3 , r\ = fi = 10 -3 , Ma = 
(3 = 1.0, kdi >l,L x = L y = L !S = 2tt. 
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Fig. 3. Inertial range turbulent spectra in kdi < 1 regime. Shown 
are the magnetic and velocity field fluctuations along with that of 
the density fluctuation. These fluctuations closely follow a k~ 5 ^'' 1 
omnidirectional spectrum in the Alfvenic regime. Simulation pa- 
rameters are the same as in Fig. (O, except di = 0.05. 



ulations [see Fig. (5)1. Applyin g the Kolmogorov' s phe 



nomenology iKolmogorov I d 1 95 lb : IKraichnan I d 1 9651) that 
the energy cascades in the inertial range are local and that 
they depend on the Fourier modes and the energy dissipation 
rate, we obtain k~ x B\ ~ [B\k 2 ) a k^ . Equating the indices 
of the common bases, we get a — 2/3 and (3 = —7/3. This 
results in an energy spectrum oc fc~ 7 / 3 that is consistent 
with our 3D simulations [see Fig. ©I and the observations. 
On the other hand, the use of tmhd hi estimating the energy 
dissipation rates retrieves the inertial range MHD spectrum. 
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This indicates that the Hall effects may be responsible for 
the spectral steepening in the solar wind plasma fluctuations 
in the kdi > 1 regime. It is worth noted that our simula- 
tions in kdi < 1 exhibit MHD like k~ 5 ^ 3 omnidirectional 



Kolmogorov-like lKolmogorov I d 1 95 It) spectrum as shown in 
Fig. ®. 



5 Conclusions 

In summary, we have presented a survey of the spectral prop- 
erties of low-frequency electromagnetic turbulence in uni- 
form magnetized plasmas. Specifically, we have focused 
on the electron whistlers, and high-frequency kinetic Alfven 
waves. While for the electron whistlers the motion of ions is 
unimportant, the ion dynamics plays an essential role in the 
high-frequency KAW turbulence. We have developed codes 
to numerically solve the governing nonlinear equations for 
the whistler mode and KAW turbulences. Our simulations 
for the whistler mode turbulence depict a dual cascade pro- 
cess between the energy and mean magnetic potential in con- 
sistent with a Kolmogorov phenomenology. Furthermore, 
the high-frequency kinetic Alfven regime, typically observed 
in the solar wind plasma, is studied by means of 3D simu- 
lations. Our 3D simulation results exhibit that the charac- 
teristic turbulent interactions in the high-frequency regime 
evolve typically on the Hall time-scales and modify the iner- 
tial range spectra that eventually leads to the observed spec- 
tral break in the solar wind. 

The support of NASA(NNG-05GH38) and NSF (ATM- 
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This work was partially supported by the Deutsche 
Forschungsgemeinschaft through the project SH 21/3-1 of 
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